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Abstract: We use a quantum-dynamical model to investigate the optical 
response of graphene under low excitation power. Ultrafast carrier relax- 
ation processes, which play an important role for understanding the optical 
response of graphene, are included phenomenologically into the model. 
We obtain analytical solutions for the linear and third-order nonlinear 
optical response of graphene, and four-wave mixing in particular. This 
theory shows agreement with recently reported experimental data on linear 
complex optical conductivity and four-wave mixing, providing evidence for 
ultrafast quantum-dephasing times of approximately 1 fs. 
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1. Introduction 

Graphene is a two-dimensional material consisting of a monolayer of carbon atom forming 
hexagonal lattices [1]. Much has been learned since its first isolation from graphite and deposi- 
tion onto a dielectric substrate [2]. Due to its linear and massless band structure near the Dirac 
points, graphene has been shown to posses astonishing linear electronic and optical properties. 
Electron mobility in graphene reaches 2 x 10 5 cm 2 /V s, much higher than normal materials such 
as silicon [3]. Graphene is also the only known material exhibiting anomalous quantum Hall 
effect under room temperature [4]. The optical absorption per layer of graphene is related to 
the fine-structure constant by na — 2.3% over a broad range of terahertz and optical wave- 
lengths [5-7]. Potential applications include ultrabroadband fast detectors, and the replacement 
of transparent conductive oxides with graphene. Also of interest are potential coherent elec- 
tronics applications that would use quantum interference of electrons to provide new device 
functionality. 



Besides linear electronic and optical properties, the third-order nonlinear response of ma- 
terials is useful for certain applications such as wavelength conversion [8] and quantum- 
entanglement generation [9]. The most important physical parameter for nonlinear-optical ex- 
periments is the nonlinear coefficient, typically the nonlinear susceptibility, of the material used. 
However, in practice, nonlinear-optical experiments are limited by other factors, in particular, 
the phase-matching conditions arising from wavelength-dependent refractive indices. The accu- 
mulation of phase-mismatch as the interacting waves propagate through the nonlinear medium 
sets a limit on the bandwidths of the nonlinear applications. Because the graphene thickness is 
much less than an optical wavelength, non-critical phase-matching occurs simply, though at the 
cost of some optical absorption. 

Nonlinear properties of graphene have been investigated theoretically in [10-14]. Theoretical 
results predict high nonlinearity of graphene at the terahertz-frequency [10, 11] and optical- 
frequency [14] ranges. More recently, the nonlinearity of graphene in the visible-optical range 
was experimentally measured by four-wave mixing [15]. When taking into account the nonlin- 
ear response per unit thickness, the theory adopted in Ref. [15] shows that graphene's effective 
third-order nonlinear susceptibility is of 8 orders higher than glass. Ref. [15] estimated 
the % °f graphene by evaluating the nonlinear surface-current density induced by an exter- 
nal electromagnetic field and estimating the dynamic nonlinear response of electrons coupled 
with an external field is characterized by the electron-photon interaction Hamiltonian. Instead 
of the dipole interaction Hamiltonian utilized by [15], a more general form of the interaction 
Hamiltonian can be obtained by the minimal substitution p — > p — qA(t), giving a more accurate 
description of the dynamics. As we will see in this paper, the dipole interaction Hamiltonian 
is not adequate to describe the quantum dynamics of electrons in graphene. Additionally, the 
quantum dynamics of electrons in graphene is also affected by ultrafast many -body interactions 
such as electron-electron and electron-phonon scattering [16, 17, 19], which in turn may depend 
on carrier densities, conduction band energies, and temperature in non-trivial ways, resulting 
in ultrafast quantum dephasing. Other quantum dephasing mechanisms may exist. A complete 
physical model needs to include these effects, which have not yet been included in previous 
theoretical calculations of the nonlinearity of graphene. 

In this paper, we present a quantum-dynamical model for investigating the quantum dy- 
namics of electrons in graphene. The electron-photon interaction Hamiltonian is obtained by 
the minimal substitution on the free-electron Hamiltonian to relate the quantum dynamics of 
electrons with the vector potential of the electromagnetic field, as in [12]. This paper differs 
in that it includes electron-electron and electron-phonon scattering by introducing two phe- 
nomenological decay rates. This paper also focuses on nonlinear mixing rather than on the lin- 
ear decay-free, dephasing-free model. In a quasi-continuous-wave pump experiment, in which 
the electron-photon coupling time (order of ps) is much longer than the carrier-relaxation time 
(order of fs), these phenomenological decay rates allow us to avoid dealing with complicated 
microscopic quantum-mechanical calculations and to obtain analytical solutions for quantum 
dynamics of electrons in graphene by exploiting quantum-perturbation method. The proposed 
model is appropriate for excitation power well below the saturation threshold of graphene (4 
GW/cm 2 ) [20], in which case the population inversion for each quantum state in the Brillioun 
zone does not significantly change. Once the quantum dynamics of electrons in graphene is 
obtained, both linear and nonlinear optical conductivity produced by different optical processes 
will be derived. 

This paper is organized as follows. In Sec. 2, we discuss the quantum dynamics of electrons 
in graphene. Our discussion in this section covers the electron-photon interaction and ultrafast 
quantum dephasing caused by carrier-relaxation processes in graphene. In Sec. 3, we derive the 
surface-current density in graphene. We calculate both the linear and nonlinear surface-current 



densities that are caused by the driving optical frequency fields. We also discuss the linear 
complex optical conductivity and connect it to ultrafast quantum dephasing in graphene. In 
Sec. 4, we apply the nonlinear surface-current density obtained in Sec. 3 to evaluate four-wave 
mixing in graphene, and fit the theory to the recently reported experimental data [15]. In Sec. 
5, we discuss the validity of the proposed quautum-dynamical model under certain physical 
circumstances. Conclusions are reviewed in Sec. 6. 

2. Quantum dynamics of electrons in graphene 

In this section, the quantum dynamics of electrons in graphene will be reviewed, starting from 
first principles. We begin with the free-electron Hamiltonian, obtained for the 2D tight binding 
model of nearest neighbour interaction approximation. The electron-photon interaction Hamil- 
tonian is derived later by minimal substitution. We then discuss the ultrafast carrier relaxation, 
which plays an important role in the quantum dynamics of electrons in graphene. The resulted 
ultrafast quantum dephasing in graphene is included by addition of phenomenological decay 
rates into the model. 

2.1. Free-electron Hamiltonian and electron-photon coupling 

In the 2D tight-binding model, or nearest neighbour interaction approximation, the motion of 
electrons is limited by assuming that they can only hop to their nearest neighbours. The Hamil- 
tonian of electrons in graphene under this assumption is written in the second-quantized lan- 
guage [1]: 

H = -V £ (al .b aJ + H.c.), (1) 

where rj w 2.8eV is the hopping energy and the product of the creation and annihilation oper- 
ators a a p a j illustrates a process in which an electron with spin a € {til} is annihilated on 
site Rj of sublattice A and a new electron of the same spin is created on its neighbour site Rj 
of sublattice B. The Fourier transform of the Hamiltonian in Eq. 1 allows the derivation of the 
free-electron Hamiltonian in the momentum space, which is 

#=L(4^)#oQ), (2) 

k 

where 

H ° = ( K o" ) (3) 

is the first-quantized free-electron Hamiltonian, characterizing the quantum dynamics of a sin- 
gle electron with momentum p = hk. The summation in Eq. 2 is performed over the entire 
Brillioun zone. In the vicinities of two Dirac points in the Brillioun zone, we have the 
following linear dispersion relation for h p : 

h P = Vf(Px T iPy) = v F p(cos 9 =F /sin 0), (4) 

where vp sa 10 6 m/s is the Fermi velocity of electrons in graphene. Energy eigenstates and 
corresponding eigen-energies can be found by solving the 2D Dirac equation [21]: 



H \(p)=E\<p). 



(5) 



Given the electron momentum p, one finds two energy eigenstates with opposite eigenvalues 
E C (y\ = ±|/j p | = ±\'fP, corresponding to the conduction-band and valence-band respectively: 



1 



'9i - \-hy\h t 

Having discussed the Hamiltonian for a single free electron in graphene, we next derive 
Hamiltonian of an electron coupled to an external electromagnetic field, which is characterized 
by its vector potential A(f ), related to the electric field E and the magnetic field B by 

. = -£ 

B = V x A. 

Suppose at time t = 0, we turn on a homogeneous a.c. electromagnetic field having a vector 
potential with two polarization components: 

A(t) =A x (t)±+A y (t)f = (X(f),A,(f)) , (9) 

where x and y are two orthogonal unit vectors of the Cartesian coordinate. Each polarization 
component A q (t), q € is the summation of different frequency modes on a polarization: 

A g (t) = £A <? ,„ e -"V' r = £A (/ ,„(f). (10) 

n n 

Both positive and negative frequencies are allowed. The complex conjugate relation n = 
A* _ n holds to guarantee a real field. By using Eq. 8, the electric field on polarization q is 
related to the vector potential by 

E q (t) = -^P- =Yj^A^~ imq ' nt = ^E qjl e-^"'. (11) 
The electric-field amplitude and the vector-potential amplitude have the following relation: 

Eq.n = i(&q,nAq,n- (12) 

The electron-photon interaction Hamiltonian can be obtained by the minimal substitution 
p — > p + eA(t ) on the free-electron Hamiltonian Ho [22] . The Taylor expansion of the electron- 
photon interaction Hamiltonian around p reads 

Due to the linear band structure of graphene, higher-order terms in the Taylor expansion vanish 
in the vicinities of two Dirac points. Let H(t) =Hq + V(t), so that we have 



VAp-eA(0 \_^r fA 



dhr 



where u e {x,y}. In the representation of energy eigenstates |V p ) and \C P ), V(t) is organized 
into the following matrix form: 



v(t) _l (Vplv(OIVp) (v P \v(t)\c p ) \ y l Vfy v$ c \ 



(C p |V(r)|Vp) (C p \V(t)\C p ) J tj '\ V cv Vcc 
Around the two Dirac points K*, by use of Eq. 4, one obtains 

Vyy = —Vf COS 9 , Vy c = ±lVf SHI 9 , Vgy = T' V F Sill 9 , Vqq = Vf COS 9 

V vv = — VF&m9,V vc = iv p cos 9 ,V^y = ±/vf cos 0, V^ c = Vfsin0. (16) 

In general, the quantum dynamics of an electron can be obtained by solving the time- 
dependent Schrodinger equation: 

ih^ = [H(t),p], (17) 
with p the density operator of the electron. 
2.2. Electron relaxation 

Eq. 17 describes the quantum dynamics of an electron with momentum p that is subject to an ex- 
ternal electromagnetic field characterized by its vector potential A(f ). A valence-band electron 
can be excited, and become a conduction-band electron by absorbing a photon. The transition 
happens with a particular probability that depends on the field intensity and photon energy. The 
excited conduction-band electron could be stimulated into the valence band and emit a photon. 
This whole process is the Rabi oscillation. However, in a complete physical picture, electrons in 
the conduction band undergo ultrafast carrier relaxation processes caused by the faster electron- 
electron scattering followed by the slower electron-phonon scattering [16-18]. These scattering 
processes destroy the quantum coherence of electrons and result in quantum dephasing. In a 
strong excitation regime, a strong pump pulse depletes the population in the valence band and 
excites electrons into the conduction band. Thereafter, the hot electrons in the conduction band 
undergo an ultrafast electron-electron scattering process (as do the holes in the valence band) 
that results in quasi-thermalized statistics for electrons in the conduction band and holes in 
the valence band. The quasi-thermalized distributions establish new state-occupation probabil- 
ities within the two bands. As a result, a strong pump not only changes the state-occupation 
probability of the state having the corresponding energy, but also affects the state-occupation 
probabilities over the whole band. A complete treatment that would include the electron-photon, 
electron-electron, and electron-phonon interactions requires sophisticated quantum-mechanical 
calculations [17, 23]. There may be additional quantum dephasing processes. However, if we 
limit ourselves to low-excitation quasi-continuous-wave regime, in which the state-occupation 
probabilities over the whole bands do not change significantly by the pump nor in time, the 
situation is much simpler. To model the ultrafast quantum dephasing in a low excitation quasi- 
continuous-wave regime, we simply introduce two phenomenological decay rates F\ and I^, 
describing the population and quantum-coherence damping that phenomenologically represent 
the overall effect of the ultrafast scattering processes. With the two phenomenological decay 
rates, one arrives at the Bloch equations for graphene under low excitation power: 



P = -Ti{p-p 



2ei 



Y,A uJ (t)[VZc(pvc + Pcv)] 



p vc = (-i(Ovc-^2)Pvc- T Y. A u.j( t )[ 2v vvPvc + yvcP] 



"■J 
ei 
h 



Pcv = (-ico C v-r 2 )pcv-^Y< A u.j(t)[-2VvvPcv + V{; c p}, 



(18) 



where p = pec — Pvv is the population inversion, fflfcy = — fiVc = 2Ec/h, and p eq is the pop- 
ulation inversion in thermal equilibrium. Compared to the Bloch equations for atom vapours 
and normal semiconductors [17, 24], both diagonal and off-diagonal driving terms exist for the 
off-diagonal elements of the density matrix. General solutions of Eq. 18 are difficult to find. 
Instead, we seek perturbative solutions, which are are good approximations when the pump 
power is well below the saturation threshold, which is measured to be 4 GW/cm 2 [20]. 

The perturbative method gives the solution in the following cascaded form in the frequency 
domain, with each order denoted by a positive index k: 



,(*+*: 



(A) 



ey 2 A 
lX. (£2-iTi 



V VC 



p [ v k] c (H + co uJ ) + p { c k j (£2 + co uJ ) 



A uJ 



2V$ v p$ (£2 + (0 lLj ) + V^cP w (£2 + co uJ 



2V$yp W (£2 + (Ouj) + V^ c p W (£2 + fflyj) 



(19) 



where £2 = co — iK with K > in order for the Laplace transform to converge. 
3. Surface-current density in graphene 

In this section, the quantum dynamics of electrons will be derived perturbatively. We first derive 
the first-order solution, which determines the linear optical response of graphene. We then 
compare our theoretical result with the linear complex conductivity measured by [25-27]. The 
cascaded relation in Eq. 19 allows the derivation of higher-order density matrix terms, which 
determine the nonlinear quantum dynamics of electrons in graphene. We continue to calculate 
up to the third-order terms, which allows the prediction of nonlinear optical responses such as 
four- wave mixing. We finally fit the theory to recently reported four- wave mixing experimental 
data [15]. 

3.1. Linear optical conductivity 

The zeroth-order solution to the density matrix is simply the thermal-equilibrium values of the 
population inversion and quantum coherence. At temperature T = 0, we have p^l = Py^. = 



Pcv = ^' Pvv = !• Inserting zeroth-order solution into Eq. 19 leads to the first-order solution: 



p (1) (£2) = 



W.lt 



where A w .„ (£2) is the Laplace transform of the time-dependent vector potential of one frequency 
mode indexed by n on polarization w. The expected velocity on q polarization of a field-coupled 
electron is given by 

Tr[v«p], (21) 



v 



where the velocity operator v q is defined as 



\+'ZMt)[ gv W "- ]=•" + /". (22) 



To simplify our discussion, we only consider an external field polarized in the x direction with 
vector potential A(t). 

The linear electron velocity is contributed by both the diagonal elements from p (°) and the 
off-diagonal elements from pW [21, 29]: 



(v v )/ = Tr 



(23) 



The linear surface-current density 7/ (£2) can be obtained by an integral over the whole Brillouin 
zone. Since the integrand (v x )i is only appreciable in the vicinity of the input optical frequency 
ft), instead of an integral limited in the Brillioun zone, we can extend the integral to infinity 
without significantly affecting the result. The integral is written as 



Ji (Q.) = -Lr / dO I -e(v x )[kdk 



[■2k 

4n 2 . 

e 2 

= (n-/r 2 )A(£2). (24) 

By using the frequency-domain relation E(£l) = — iQA(£i), the linear optical conductivity is 
derived as ^ 

a i= 4x T^( l + — l ) =Oo(l + T-i), (25) 
I6h\ (o J \ co J 

where Oq is the theoretical prediction of the universal optical conductivity of graphene [5- 
7]. A factor 4 is added due to valley and spin degeneracies. If we compare Eq. 25 with the 
experimental measurements of the complex optical conductivity of graphene at A = 550nm in 
[25-27], we find that T2 ranges from 1.39 x 10 15 s _1 to 4 x 10 15 s _1 by normalizing the real 
parts to Co- 
in [28], the authors ascribe the imaginary part of the conductivity of graphene to the vir- 
tual transition of electrons at M and Y points in the Brillioun zone. However, our quantum 
mechanical calculations based on the nearest-neighbour-interaction approximation find that we 
have V/ip = at the Y point, leading to a zero interaction Hamiltonian, i.e., V(t) = 0. At the 



M point, the interaction Hamiltonian gives pure real diagonal elements, and pure imaginary 
off-diagonal elements in Eq. 15. Thus, the conductivity at the M point is still real as long as 
we set r"i = Y2 = 0. Thus it seems that the relatively large imaginary part of the conductivity 
arises from the ultrafast quantum dephasing in graphene. The fact that this time constant is ap- 
proximately 1 fs indicates both that optical measurements can be used to characterize quantum 
dephasing in particular samples and that dipole quantum coherence would seem to be a poor 
candidate for coherent electron devices in graphene. 

3.2. Nonlinear optical conductivity 

To evaluate the nonlinear interactions in graphene, we need to solve for higher-order density- 
matrix terms p( 2 ' and p^\ Having obtained p^ and p' 3 ', the expected nonlinear velocity on 
polarization q for an electron with momentum p can be calculated as 



(v q )nl = Tr 



vV 3 ) 



(26) 



According to Eq. 22, both p( 2 ' and p^ contribute to (v ? )„/. To derive the nonlinear surface- 
current density, we need to integrate over the entire Brillouin zone. p( 2 ' and p( 3 ' both are 
appreciable only in the vicinities of their resonant frequencies. The contribution from p^> can 
be neglected since the second derivative of h p gives zero due to the linear band structure of 
graphene near the Dirac points. Thus, in the nonlinear surface-current-density calculation, we 
only keep the v' q term and drop the v" q term in Eq. 22. The resulting nonlinear velocity is 
written as 

EEE >l /;;(; , ; , ; „ . (27) 



= Tr 



vV 3) 



where the expected velocity of one frequency component on polarization q is defined as 
(v q , ) . Since the main contribution of (v q , ) comes from the vicinities of each res- 
onant frequency, we can extend the integration in the Brillouin zone to infinity and without 
significantly affecting the result. By assuming the Dirac cone goes to infinity, we obtain analyt- 
ical solutions to the third-order surface-current density of at one frequency component having 
polarization q to be 



2k 



= -l-^AuiA 



16ft 



^ v M,/;v,m;w,/i 



)kdk 



■ + 



(£2 + <»„,/ + -ir 2 ) 



£2 - ;Ti (£2 + 0) u j — 1T1) (2£2 + co u j + ffl,,,,, - 2/r 2 



(£2 + ©„,/ + ©,,„,), 



(28) 



where M?,«,v,w and V q ,u.v.w are 



nvp Jo 



2k 



V$ v V$ c VZyVy c d6 



2k 



Vq,u,v,w — —4 / VvC V VC V VC V Vcd® 



(29) 



In the time domain, the steady-state solution to the nonlinear surface-current density reads 



Vp€ E u [E vm E wn 

16fj 3 ft)„,/ft),.,,„ft),y,n 



Mc/.M.l 



((0 uJ + ok,, 



V 9 ,u,v,w(fl\v,B + lT2) 



(C0u,l + «v,m + 2<0 Ky) + 2;T 2 ) (ft),.,,,, + ft),,.,,, + (Ti ) 



■CDwst + iTi) 

r,-i('"iiJ + 0>v,m + 0} w ,n)t 



(30) 



where a coefficient 4 is added to include valley and spin degeneracies. We have substituted the 
vector-potential amplitude with the electric-field amplitude by use of Eq. 12. Having obtained 
the surface-current density of one frequency component/^, (f ) on polarization q, the total 
surface-current density on this polarization is the summation over all frequency components, 



^(0=LEE^ ; , w (0. 



(31) 



u.l v,m w,n 



Eq. 3 1 is a general expression for the total third-order nonlinear surface-current density, com- 
posed of different frequency modes, each of which is produced by a particular nonlinear-optical 
process. For example, let co Pl , co P2 be the frequencies of the pump modes and co s be the fre- 
quency of the signal mode, the frequency 3ft) Pl results from third-harmonic generation of the 
pump at frequency (0 Pj while ft), = ft) Pl + ft) p , — ft), is produced by four-wave mixing. 

We mark one frequency mode with frequency ft) on polarization q as J%(f). 4(0 is the 
summation over all frequency components J q u i. vnrwn {t) satisfying ft),,,/ + ft),.,,,, + ft),,,,, = ft): 



4(0 



E 



f (t) 

u,l;v,m;w,n V ' 



Tip' 



(32) 



4. Four-wave mixing in graphene 

Of particular interest for this paper is the study of four-wave mixing, a nonlinear-optical process 
involving four modes. Two of these serve as the pump modes, one as the signal mode, and one 
as the idler mode. The optical frequencies co Pl and co P2 correspond to the two pump modes, co s 
to the signal mode, and ft), to the idler mode. Energy conservation gives co Pl + (0 P2 = ft)., + ft),. 
Microscopically, four-wave mixing is a process in which two pump photons are annihilated 
while one signal and one idler photon are created. In the special case where the two pump 
modes share the same optical frequency, i.e., co pl = co pj = CO p , the four-wave mixing is pump 
degenerate. Since four- wave mixing is a third-order nonlinear process, its strength in the electric 
field is proportional to the nonlinear susceptibility j( 3 '. If the frequency of the optical modes 
are far away from the resonant frequency of the nonlinear medium, the results is a non-resonant 

(3) (3) 

susceptibility Xnr - ^ n tne quantum-mechanical picture, %nr results from the fact that the photon 
energies are far from the bandgap of the nonlinear medium. A typical value for the non-resonant 
susceptibility in optical fibers is Xnr ~ 10~ 15 esu [24]. However, a much higher %^ [30, 31] 
is present if the bandgap of the nonlinear medium is close to the interacting-photon energies. 
The resonance-enhanced nonlinear susceptibility Xr^ is typically of the order of 10~ 7 esu. 

Graphene is a zero-bandgap semiconductor with linear band structure near the Dirac points, 
making it "more resonant" than typical resonant media. The zero bandgap of graphene results 
in the fact that any frequency within the range from DC to optical frequencies is resonant. This 
is the physical origin of the high and uniform absorption of graphene. While absorption is a 
resonance-enhanced linear process, nonlinear processes can also be enhanced by resonance. 
In particular, four-wave mixing in graphene can be enhanced by the five resonance-enhanced 
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Fig. 1. Five resonance-enhanced four- wave mixing processes in graphene. (a) Two-photon 
absorption enhanced four-wave mixing, (b) One-photon absorption enhanced four-wave 
mixing, (c) Idler resonance-enhanced four-wave mixing, (d) Signal resonance-enhanced 
four-wave mixing, (e) Detuning-enhanced four-wave mixing. 



processes plotted in Fig. 1. Different resonance-enhanced four- wave mixing processes exist 
for electrons with different energies in graphene. We next discuss these resonance-enhanced 
four-wave mixing processes. 

(a) Two-photon absorption enhanced four-wave mixing: Two-photon absorption happens 
when the bandgap of a material is twice the photon energy of the incident light. In graphene, 
this corresponds to (Ocv —2co p . 

(b) One-photon pump absorption enhanced four-wave mixing: Unlike two-photon absorp- 
tion, in which one electron is excited by two incoming photons simultaneously, an electron is 
excited by a single photon in a one-photon absorption process. In graphene, the condition for 
one-photon absorption is (Ocv — (0 p . 

(c) Idler-enhanced four-wave mixing: In four-wave mixing, if the emitted light frequency, 
i.e., (Oj, is close the the bandgap of the nonlinear medium, we also have resonance-enhanced 
effect. In graphene, this happens at the bandgap a>cv = GH- 

(d) Signal-enhanced four-wave mixing: If the signal-photon energy is close to the bandgap 
of the nonlinear medium, we also have enhanced effect for four-wave mixing. In graphene, 
signal-enhanced four-wave mixing happens at bandgap (Ocv = ®s- 

(e) Detuning-enhanced four-wave mixing: Detuning-enhanced four-wave mixing in 
graphene is similar to the four-wave mixing enhancement by coherent anti-Stokes Raman 
scattering [32, 33], when the frequency detuning between the pump mode and the Stokes mode 
is close to the molecule-vibration energy. In graphene, this is the case when the bandgap is 
close to the frequency detuning between the pump and the signal, i.e., <x>cv — ©p — ©s- 



4.1. Non-degenerate four-wave mixing 

AAA. Co-polarized non-degenerate four wave mixing 

In the co-polarized case, we let all frequency modes be x polarized. We define the conductivity 
of the co-polarized non-degenerate four-wave mixing to be 

% N (t) = a C NE Pl E P2 E:.e- im -'+c.c. = /£,.(/) (33) 

The frequency components contributing to Gcn are listed in Tab. 1 . We insert terms in Tab. 1 



Table 1 . Frequency components for co-polarized non-degenerate four- wave mixing 
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into Eq. 30 and Eq. 32 to obtain 



Ocn = F "I"' p » ■ (34) 



4.1.2. Cross-polarized non-degenerate four wave mixing 

For cross-polarized non-degenerate four-wave mixing, two pumps modes at frequencies co lh 
and co Pv are on orthogonal polarizations x and y. We assume that the signal mode at frequency 
(0 S is polarized on x. Momentum conservation guarantees that the idler mode is produced on 
the y polarization. We define the conductivity for this nonlinear process to be 

J y XN {t) = o XN E Px E Py E*e- iai < +c.c.=J y eH {t) +J y _ a ,(t). (35) 

The contributed frequency components to Gxn are listed in Tab. 2. By comparing Table 2 with 



Table 2. Frequency components for cross-polarized non-degenerate four-wave mixing 
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Table 1, we find that Gxn = \ &cn- Thus, Gxn and Gcn differ only by a factor of 3 as in typical 
nonlinear media. We can adopt the same approach as in [15] to compare the effective third-order 
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Fig. 2. The ratio between the third-order susceptibility of graphene and and that of glass. 

(3) 

Xgr is calculated from the nonlinear optical conductivity of non-degenerate co-polarized 
four- wave mixing for different pump detunings and decay rates Ti and Ti- Pq denotes the 
center wavelength between the two pumps. 



susceptibility of graphene with the third-order susceptibility of silica. The effective third-order 
susceptibility of graphene is given in the MKS units by 

^ = ^4, 06) 

where d ~ 3. 3 A is the effective thickness of a monolayer of graphene. Typical value for non- 
resonant third-order susceptibility of silica such as glass is %ty = 2.22 x 10~ 23 m 2 /V 2 [24]. We 
plot the ratio between the two susceptibilities in Fig. 2. In each figure, three curves correspond 
to center pump wavelengths at 550nm, 775nm and 1550nm respectively. For each center pump 
wavelength, we choose three detuning frequencies at lOTHz, 20THz and 30THz between the 
two pump modes. 

Although Ref. [15] predicts that the nonlinear susceptibility of graphene is approximately 8 
orders of magnitude greater than that of insulators, our more complete model shows that the 
ratio between the two susceptibilities can vary from 5 to 9 orders of magnitude, depending 
on the pump wavelengths, the signal detuning, and the decay rates T\ and Ti- A faster quan- 
tum dephasing results in a lower nonlinearity due to the fact that the coherence of electrons is 




Fig. 3. The four-wave-mixing current intensity ratio between Eq. 37 and the theoretical 
result in [15]. Five curves corresponds to five different pump and signal wavelengths in 
[15]. We set F\ = 2F2- Both the ratio and the decay constants are plotted in logarithm 
scale. 



quickly destroyed by the environment via scattering processes. Our theory also predicts that ss 
the pump-signal detuning becomes smaller, the strength of four-wave mixing increases due to 
resonance-enhanced effect. Another interesting phenomenon is that as the quantum dephasing 
times approach 1 fs, the typical symmetry in the nonlinear susceptibility breaks down. 

4.2. Degenerate four-wave mixing 

The optical conductivity for degenerate four-wave mixing can be obtained by setting co Pl = 
(0 P2 = C0 P in the derivation of the non-degenerate four-wave mixing conductivity and dividing 
the result by 2 due to degeneracy. In experiment, the observed idler photon-current intensity / 
relates to the degenerate four-wave mixing conductivity Go by 

/-|£;£> D | 2 , (37) 

with Op as a function of (0 P , C0 S , F\ and Fx- In Fig. 3 we compare the idler photon-current 
intensity predicted by Eq. 37 with that is predicted by Eq. 4 in [15]. Depending on the de- 
cay rates F[ and T^, the ratio between the two current intensities varies. The result shows that 
faster quantum-dephasing rates lead to a weaker four-wave mixing. For decay rates at 10 15 
s _1 , the four- wave-mixing photon-current intensity is roughly equal to what was predicted in 
[15]. We next fit our theory to experimental data obtained by [15]. We use three values for F2 
found in Sec. 3. 1 by investigating the imaginary part of the linear complex optical conductivity 
measured by experiments. Thus, the only parameter to determine is F\. Fits to experimental 
four-wave-mixing data yields Fi ranging from 0.8 x 10 15 s to 1.25 x 10 15 s _1 . The fits are 
shown in Fig. 4. To summarize, the decoherence time is taken from linear refractive index 
measurement, without fitting. Because the amplitude of the nonlinear response is taken directly 
from the theory at the decoherence time, only one parameter is used to make the fit. The two 
decay rates were both obtained from experimental data, which provides evidence for ultrafast 
quantum dephasing in graphene at the 1 fs timescale. Our fully quantum-mechanical model 
results in faster decay rates than what was reported in [20] based on semi-classical calculations. 
We note here that the decay rates we found are different than carrier-relaxation time [16, 20] 
in two ways. First, carrier-relaxation time was obtained by differential transmission (DT) ex- 
periments where the graphene was excited by an ultrafast strong pump and then probed by a 
weak pulse. In the probing phase, only electron-electron and electron-phonon interactions ex- 
ist. However, in a continuous-wave or quasi-continuous-wave experiment, the relatively strong 
pump (still weaker than the saturation threshold) interact with the graphene through the whole 
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Fig. 4. Fitting of the theory to the experimental data obtained by [15]. The two phe- 
nomenological decay rates for the three curves are Ti = 0.8 x 10 15 s~ 1 , TS = 4 x 10 15 s~ 1 , 
Ti = 0.9 x 10 15 s~\ T 2 = 2.69 x 10 15 s -1 , and Ti = 1.25 x 10 15 s -1 , T 2 = 1.39 x 10 15 



process. Thus, we expect the decay rates to be faster since electron-photon, electron-electron 
and electron-phonon interactions all happen simultaneously. Second, DT experiments measure 
the population inversion at the given time. However, our decay rates also include the damping of 
quantum coherence of electrons. In addition, pure dephasing processes may exist in graphene. 

5. Discussions 

Finding a correct Hamiltonian for a quantum system is of course the key step in derive its 
dynamics. As we have discussed, we derive the Hamiltonian of electrons in graphene by the 
minimal substitution p — > p qA. The Hamiltonian obtained by minimal substitution is gen- 
eral and it produces the correct electron-photon interactions. However, sometimes the dipole 
interaction Hamiltonian/// = er E was used for semiconductors [17] and atom vapours [30]. 
The dipole interaction Hamiltonian can be derived from the original electron Hamiltonian by 
the minimal substitution to get 

H=^t + V(r), (38) 
followed by replacement of p by means of the commutation relation [17] 

[r,Ho] = —> (39) 
m 

We note here that Eq. 39 only holds for parabolic energy-dispersion relations. However, the en- 
ergy dispersion is linear near the Dirac points in graphene. Thus, the dipole interaction Hamilto- 
nian does not apply to graphene in the area near the Dirac points. Near the M point and F point 
in graphene, we do have parabolic energy-dispersion relations. If we compare the graphene 
Hamiltonian in Eq. 13 with the typical Hamiltonian in Eq. 38 for normal materials, we find 
that the Hamiltonian of graphene commutes with the momentum operator, and therefore we 
have energy-momentum degeneracy in graphene. At the M point and the T point, the transition 
matrix element (Cp|p|V p ) = 0. Thus, as previously discussed, electrons near the F point do not 
interact with the field in the tight-binding model within the nearest neighbour interaction ap- 
proximation. Electrons near the M point only contribute to the real part of the conductivity if 
the quantum-dephasing is not present. 

The perturbative approach is a good approximation when the saturation effect can be ne- 
glected. When saturation effect starts to become significant, the population inversion p be- 
comes so large that our perturbative solution can fail [34]. However, as we have shown, there 



exist five different resonance-enhanced four-wave mixing processes in graphene. We find that 
the population inversion p^ approaches to its maximal value around co C v — (Op, induced by 
the one-photon absorption. Thus, the one-photon absorption enhanced four-wave mixing is 
the most likely to saturate. If we inspect the purturbative solution p( 2 \ we find that popula- 
tion inversion induced by one-photon absorption is a second-order process. In comparison, the 
population inversions induced by two-photon absorption, pump-signal detuning, or idler simu- 
lation are third-order processes with much weaker magnitudes. Although the signal mode can 
also induce second-order population inversion, it can be neglected when the strong pump ap- 
proximation is valid, since the signal intensity is much less than the pump intensity. Under high 
excitation power levels that can completely change the state-occupation probabilities over the 
valence and the conduction band, a full-quantum-mechanical many-body theory is required to 
derive the quantum dynamics of graphene [17, 23]. The quantum many-body theory can also be 
used to calculate the instantaneous response of electrons in graphene. For steady-state electrons 
under high excitation power, a quasi-thermalized carrier statistics is established within both the 
valence and the conduction band. Knowing this quasi-thermalized carrier statistics, one could 
plug this into p eq in Eq. 18 and derive the steady-state quantum-dynamics of electrons under 
high excitation power. 

In our calculations, we made a linear approximation of the band structure of graphene near 
the Dirac points. This approximation is valid up to 4eV of bandgap given the hopping energy 
77 ps 3 eV. For photon wavelengths below 500nm, the approximated linear band structure starts 
to be inaccurate and the optical conductivity diverges from the universal value Gq. A full-band 
calculation [35] would be needed for this situation. 

6. Conclusions 

We have used a quantum-mechanical model to investigate the quantum dynamics of electrons in 
graphene. The tight-biding model with nearest neighbour approximation and linear band struc- 
ture near the Dirac points are assumed in order to obtain analytical results. We have also in- 
cluded ultrafast decay and quantum-dephasing effects by adding two phenomenological decay 
constants into the model. Matches to linear complex optical conductivity and fits to nonlinear 
four-wave mixing experiments establish quantum-dephasing rates of approximately one fem- 
tosecond. This provides evidence for ultrafast dipole quantum-dephasing in graphene. Our ap- 
proach shows that continuous wave or quasi-continuous-wave optical measurements can yield 
information about ultrafast scattering processes in graphene. This work also has relevance for 
proposed coherent electronics applications. 
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